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A semiclassical description of structure functions in DIS at small x is 
presented. It gives an intuitive picture of the transition from the Dou- 
ble Leading Log approximation at large Q 2 , to the power like dependence 
on x in the BFKL region at limited Q 2 . Formal derivations from pcrtur- 
bative QCD, e.g. in the BFKL or the CCFM formalisms, are technically 
complicated, and therefore such an intuitive picture may be valuable and 
possibly helpful in the work towards a better understanding of DIS, and of 
the strong interaction in general. 

PACS numbers: 12.38-t, 13.60.Hb 

1. Introduction 

Experiments which study deep inelastic ep scattering (DIS) have given 
decisive information about the basic structure of matter. Results from the 
linear accelerator at SLAC [Tj gave the first evidence for a pointlike substruc- 
ture in the proton. At the higher energies available at the HERA collider 
at DESY it has been possible to penetrate still deeper into the proton. DIS 
has the advantage over e.g. e + e~-annihilation and hadronic collisions, that 
there are two different scales, the total hadronic energy W and the photon 
virtuality Q 2 , which can be varied independently. This offers a unique possi- 
bility to study the interplay between the perturbative and non-perturbative 
regimes in strong interaction. 

A remarkable result at HERA is the observation that for a fixed limited 
value of Q 2 , the 7*p cross section grows very rapidly for large W or small 
XBj = Q 2 /(W 2 + Q 2 ), approximately proportional to x~j^ ~ (W 2 ) x with 
A ~ 0.3. Such a strong increase, which is significantly stronger than the cor- 
responding growth of the pp cross section at similar energies, was predicted 
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by the BFKL [2] evolution formalism, based on the exchange of perturba- 
tive gluon chains. Indeed, in leading logarithmic order the BFKL equation 
predicts an even stronger increase than what is experimentally observed. 

The significance of the BFKL result is somewhat reduced by the facts 
that first the next to leading order correction is very large, secondly the 
introduction of a running coupling would favour chains close to the non- 
perturbative region, which makes the result sensitive to a necessary soft 
cutoff. Both of these effects make the perturbative expansion less reliable. 
Although the experimental results do indicate that some hard perturbative 
mechanism is at work, it has also been possible to fit the data without the 
BFKL mechanism, provided the soft nonperturbative input gluon distribu- 
tion grows sufficiently rapidly for small x Bj • Measurements of the total cross 
section, described by the structure function F2, cannot separate between the 
different mechanisms. Further experimental studies of the properties of the 
final states are needed, in parallel with theoretical work to distinguish be- 
tween different possibilities. 

The derivation of the BFKL equation, and also of the CCFM equation [3] 
which interpolates between the finite-Q 2 -BFKL and the large-Q 2 regimes, 
are technically quite complicated. A semiclassical description, which can 
give an intuitive picture of the results, may therefore be valuable and pos- 
sibly helpful in the work towards a better understanding of DIS, and of 
the strong interaction in general. An attempt for such a picture will be 
presented in the sections 2-5 below. 

The Linked Dipole Chain model (LDC) [3] is a reformulation and gener- 
alization of the CCFM model. The formalism in the LDC model has great 
similarities with the description in the semiclassical model presented here, 
and thus it offers a quantitative motivation for the qualitative results in our 
intuitive picture. A very brief description of the LDC model is presented 
in section 6, including also a discussion of possible applications to hadronic 
collisions. 

In this article I want to emphasize the similarities between QCD and 
QED, and have therefore repeated wellknown results from electrodynam- 
ics. Large parts may therefore appear rather trivial to many readers, who 
consequently may omit these sections. 

2. Bremsstrahlung in Electrodynamics 

2.1. Photon emission 

Let us study a charged particle with momentum p\ , which emits a photon 
with momentum q, and after recoil has obtained momentum pi = pi — 
q, see fig. ^ The interaction between a charged scalar particle and the 



Fig. 1. A particle with momentum pi emits a photon with momentum q, and 
obtains after recoil momentum p2 = pi — q. 




Fig. 2. e + e -annihilation into a pair and a photon. 



electromagnetic field is described by the interaction term 

^ej^Aft ~ e(pi + p 2 W- (1) 

This contribution depends only on the trajectory of the charged particle, and 
corresponds to what is called the "eikonal current" . The resulting emission 
spectrum can be presented in different equivalent forms 

d6 2 dq dq 2 ± dz dq\ 

dn ~ a— ~ a — ^ ~ a — ^-dy, (2) 

9 2 q q± z qj_ 

where z = <7o/Pi,o is the fraction of the parent energy given to the photon, 
and y = iln22±k. 

a 2 qo-QL 

As an example we look at e + e~-annihilation into a pair and a 

photon, see fig. [2J If the scaled energies of the fi~ , the fi + , and the photon 
are denoted x±, x 2 , and x$ respectively, the cross section is given by 

dn ^ a (l- Xl )(l-x 2 ) ^ +X2) - (3) 
With a suitable definition of a longitudinal direction we have 

q\ = s(l - xi)(l - x 2 ), y = \ H 1 —^). (4) 

I 1 — X2 
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With this definition the factor dx\dxs/(l — Xi)(l — x$) in eq. (j3J) is identical 
to the expression in eq. ©, (dq\/q\)dy. The last factor in eq. (j2J) is a 
correction factor due to the muon spin. 

For collinear emission q± — > 0. The squared propagator gives a factor 
~ l/q±, but for such emissions we have (pi + P2) //qle, which implies a 
suppression such that the net result is proportional to the factor dq\jq\ in 
eq. ((2j). Thus we note that the vector nature of the photon is essential. 



Sudakov form factors 

The total emission probability in eq. (J3J) diverges, which thus gives an 
infinite cross section. To order a this infinity is compensated by virtual 
corrections to the probability for no emission, which is a negative contribu- 
tion. Thus in leading order a cutoff is needed in q± or in invariant masses 
M„+ 7 and M JfJL -. With such a cutoff the cross sections for zero and one 
emission are both positive, and their sum, the first order result for the total 
cross section, is approximately the same as the zero's order expression. This 
result can be generalized to higher orders. In higher orders the probability 
to emit an extra photon is also compensated by a reduced probability for 
no emission, in such a way that the total cross section is (approximately) 
unaffected. Thus the annihilation process can be separated into two inde- 
pendent phases, the initial production of the fJ. + ^~ pair, and the subsequent 
photon bremsstrahlung. 

Soft emissions, for which the recoils can be neglected, are also fully 
uncorrelated, and the probability to emit a fixed number of photons in a 
given phase space region is given by a Poissonian distribution. This implies 
that the probability for no photon emission in a certain excluded region 
exponentiates, and is proportional to the "Sudakov form factor" 



This result is quite general. For processes in which an emission does not 
change the total reaction probability, the emission probability is normally 
associated with such a Sudakov form factor. We will in the following see 
sevaral examples of this. 

Formation time 

Assume that the charged particle moves along a curved and twisted 
trajectory. The emission probability is proportional to 



2.2. Ordered cascades 




excl. region 



dn). 



(5) 




(6) 
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From this expression we see that only short wavelengths can be sensitive to 
small deviations in the trajectory on a small timescale. For long wavelengths 
the integral gives the average over a larger region in space and time. This 
averaging implies that it takes a finite time to determine that a photon has 
been emitted, which can be formulated in terms of a (Landau-Pomeranchuk) 
formation time 

r ~ l/q ± . (7) 

Ordered emissions 

The result in eq. (J2J) implies that in e.g. an e + e~-annihilation event the 
photons with large q± are emitted rapidly, directly after the annihilation, 
whereas those with smaller q± are emitted afterwards over an extended 
time period. When many photons are emitted, the ordering in q±_ also 
corresponds to an ordering in time. 

Looking for the first emission in time means looking for the emission of 
the photon with largest q±. The probability that no photon is emitted with 
transverse momentum larger than q± is determined by a Sudakov form fac- 
tor. Thus the probability for emitting a photon with transverse momentum 
q± as the first emission is given by 

dx i dx 2 / 2 2 \ n / \ f \ 

where the integral in S(q±) extends over all transverse momenta larger than 

S(q ± ) = exp(- / dn). (9) 
Jq' ± >q± 

Similarly, if first a photon with transverse momentum q± i is emitted, 
a second emission of a softer photon with transverse momentum q± 2 is 
associated with the form factor S(q± i i,q±p), defined by 

S(qi_,i,q±,2) = exp(- / dq±- — ). (10) 
Jq±,2 dq± 

This process can be repeated to give a (/^-ordered cascade of emitted pho- 
tons. 

Spacelike cascades 

In spacelike cascades we have bremsstrahlung both as initial and as final 
state radiation, as indicated in fig. |21 Final state radiation is very similar 
to the emission in e + e~-annihilation discussed above. The ordering due to 
formation time implies that also the initial radiation cascade is gj_-ordered: 

Q±,i < Q±,2 < ■ ■ ■ < q±, n < Q 2 (11) 



G 




Fig. 3. A spacclike photon cascade. 




q, £3 



Fig. 4. Triple gluon vertex. 

As mentioned earlier, neglecting recoils all photon emissions are uncorre- 
cted, and therefore the ordering between them is not relevant. The ordering 
is only important when the recoils cannot be neglected, and therefore the 
current is modified for the subsequent emissions. However, as we will dis- 
cuss more in the following, the ordering is very essential in QCD, where the 
emitted radiation is not neutral, but carries colour charge. This implies that 
an emitted gluon changes radically the current, which determines later (i.e. 
softer) radiation. Therefore the order of the emissions becomes important, 
and is an essential feature of the evolution processes. 

We also note that the separation between initial and final state radiation 
in fig. |S] is gauge dependent. This separation is not determined by nature, 
but depends on the formalism used. 

3. Bremsstrahlung in QCD, 
e+e - -annihilation 

3.1. Gluon emission 

In QCD the gluons carry colour charge, and therefore radiate, or split 
into two gluons. The dominant interaction is the triple-gluon coupling (for 
notation cf. fig.0J £3 denotes the polarization vector for the "emitted" gluon 



(a) 



(b) 



Fig. 5. A quark and an antiquark emit gluons coherently by diagrams (a) and 
(b). In the angular region between the quark and the gluon the red and the 
antired charges radiate softer gluons approximately independently, forming a rf 
colour dipolc. At larger angles the emission from the red and the antired interfere 
destructively, and the quark-gluon system radiates as a single blue charge. Thus 
in this region the emission corresponds to a bb dipolc. 

with momentum q): 

g{[(pi +P2)e 3 ] (eie 2 ) - [(q+Pi)e 2 ] (e 3 ei) - [(p 2 - q)ei] (^3)} (12) 

The first term is an eikonal contribution similar to the bremsstrahlung from 
a quark or the electromagnetic radiation from an electron in eq. ©■ The 
factor e\C2 implies that the recoiling parent keeps its polarization, and it 
gives an emission density proportional to 



This term dominates for soft emissions. The second and third terms in 
eq. H12[) correspond to non-eikonal currents, in which the polarization vector 
for the emitted gluon, €3, is multiplied by the polarization vector for the 
parent gluon, either before or after the recoil. 



Study as an example the reaction in ng.Efci, where a blue-antiblue quark- 
antiquark pair is produced, and the blue quark emits a blue-antired gluon 
and becomes red. An identical final state can be obtained from the diagram 
in fig.Et>, where the gluon is emitted from an initially red antiquark, which 
becomes blue. The two diagrams interfere, and it is not possible to tell who 
is the parent to the emitted gluon. It therefore makes sense to regard the 
quark and the antiquark as a colour dipole, which emits gluons coherently. 




(13) 



3.2. Colour coherence 



Fig. 6. A colour dipole cascade. 



The emission probability is given by the corresponding expression for photon 
emission in eq. ©, with a replaced by a s . 

In the region between the quark and the gluon, the separation of red 
and antired charges will give radiation of softer gluons. However, in direc- 
tions further away from these particles the emission from red and antired 
interferes destructively. Thus in these directions the emission corresponds 
to the blue and antiblue charges of the initial qq pair, i.e. the emission 
from a bb colour dipole. This interference effect is generally called angular 
ordering In the restframe of the quark and gluon, the emission of softer 
gluons corresponds just to a rf dipole. Thus the emission of softer gluons 
corresponds to two independent colour dipoles. 



3.3. Gluon cascades 

It is possible to describe a parton state in two equivalent (dual) ways 
[5], either in terms of momenta, and spins, Si, for the gluons, or in terms 
of momenta, pi, and directions, di, for the dipoles. Thus gluon emission 
can be described as a process in which one dipole is split into two dipoles, 
which are split into four, etc., as shown in fig. El This formulation is used 
in the Dipole Cascade Model |)J[7|. It is the basis for the Ariadne MC jSJ, 
which has been very successfully applied to e + e~ annihilation data, cf. e.g. 
ref. JUJ. 

The dipole formalism is also convenient for analytic studies of the prop- 
erties of QCD cascades, and we also note that the resulting chain of colour 
dipoles gives a very natural connection to the Lund String Fragmentation 
model [TO]. 
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4. DIS, DGLAP evolution 

4-1. Notations 

The density of quarks and gluons within the proton is described by 
"parton distribution functions" f q (x) and f g (x), where x is the fraction of 
the proton momentum carried by the parton. We will use capital letters to 
denote the densities in the variable lnl/x, and for simplicity we will write 
F(x) and G(x) to denote x- f q {x) and x-f g (x) respectively. 

The electron-proton cross section is expressed in terms of structure func- 
tions F\ and F2. In the parton model these functions are related to each 
other, and to the quark density. The exchanged photon is absorbed by a 
quark, which after the absorption stays on the mass shell. If g 7 and p p denote 
the momenta of the photon and the proton, respectively, energy-momentum 
conservation is implemented by a delta-function 5(x — Q 2 /2p p -q~ ( ), where 
Q 2 = —q 2 , and the ep scattering cross section can (neglecting less important 
kinematic factors) be written 

- a V^) = ^E^(^ (14) 



dQ 2 d\nl/x Q 4 v 1 Q 4 



<i 



where x = XBj = Q 2 /2p p -q- ( , the sum runs over quark species q, and e q is 
the corresponding quark electric charge (in units of the elementary charge). 
The quantity can also be interpreted as the total "y*P cross section. 

In the distribution functions F(x) and G(x) the parton densities are 
integrated over transverse momenta. In the following we will also discuss 
"non-integrated" parton distributions, which depend also on k±, and which 
we will denote by curly letters, JT(x,/c^) and Q(x,k 2 _) for the quark and 
gluon densities respectively. 



4-2. Ordered ladders 

In the parton model the structure function F2 and the parton distribu- 
tions are functions of x only, but in QCD they also receive a dependence 
on Q 2 . Assume that a quark with momentum pq = x^p emits a gluon with 
momentum (1 — z)po before interacting with the photon, as shown in fig. d 
This implies that x = Q 2 /2pg 7 = zxq. The contribution to the cross section 
from the eikonal current is then determined by the emission probability 

Aa s dz dq\ 

Frob. ~ +o{x — zxq). (15) 

6TT 1 — z qj_ 

When many gluons are emitted as in fig. EH we again obtain an ordered 
emission determined by the formation times: 

qli < q 2 ±2 < ■ ■ ■ q± n < Q 2 (16) 
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Fig. 7. A quark with momentum po = x aP emits a gluon with momentum (1 — z)po 
before interacting with the photon. 

Thus the contribution to the structure function F from a chain with n links 
is given by a product of n factors 



4a, dzi 



dq 2 



±i 



3tt 1 



(17) 



where x = (J[ n Z{)xq, and the transverse momenta satisfy the ordering con- 
dition in eq. ((TBI) . 

Although many gluons are emitted in the process in fig. [SJi, there is 
still only one quark, which can be hit by the photon. When Q 2 is large, 
the many gluons, which can be emitted before the jq interaction, imply a 
reduced probability to collide without gluon emission. The factor in eq. (|17|). 
which determines the probability to emit a gluon, must then be multiplied by 
a Sudakov form factor, S, which describes the probability that no emission 
has occurred between q±i—% and q±i- Thus we have to make the following 
replacement in eq. (|17|) : 



1 



9(1 



dq 



Zi q 2 ± . 



1±,i 



where 



S(q 



2 a 2 



exp 



Aa s dq 



3vr qj_ Jo 1 



dz 



(18) 



(19) 



We have here introduced a cutoff e, which in a MC must be kept different 
from zero, since e — > implies that S" — >• 0. A small e gives a more accurate 
result, but also a slower program. 



4-3. DGLAP evolution for large x 

Summarizing the results in the previous subsection, and summing over 
the number of links, n, in the chain, we find the following expression for the 
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Fig. 8. (a) A chain formed by initial state radiation, (b) To describe properties 
of exclusive final states final state radiation has to be added. Note that a soft 
final state gluon also can be emitted from a virtual link, if the recoil is small. 
Such emissions, like the one emitted from the link k-i in this figure, will be further 
discussed in section 5.5. 



structure function: 



En I / ^i^ff 5 ^'^- 1 ^^-^-^ 



n i 



x5(x-Y[z r x )6(Q 2 -ql n ) (20) 

j 

We note that in this result most of the z- values are close to 1. We also 
note that, due to transverse momentum conservation and the ordering in 
eq. (|16|) . the transverse momentum of the (quasi)real emissions, q±i, are 
approximately equal to the transverse momentum of the virtual links, which 
in fig. |SK are denoted Thus the variables q±^ in eq. 1)20(1 can to the 

given accuracy be replaced by the link variables k_\_,i- 

It is also possible to include the effects of the spin of the quarks by the 
exchange 

_L &±£L (2D 

l-z 1-z v ; 

Here the numerator corresponds to (half) the factor (a;f + x|) in eq. (|HJl. and 
does not change the behaviour near the dominating singularities at Z{ = 1. 

Taking the derivative of eq. (|20|) with respect to In Q 2 gives the following 
differential- integral equation (the DGLAP evolution equation): 



dF(x,Q 2 ) _ 4a 
dlnQ 2 



— [ dzdx'P(z)F(x',Q 2 )5{x-zx'). (22) 

3-7T J 



It has here been possible to take the limit e — > 0, and thus remove the 
dependence on the cutoff. The "splitting function" P in eq. (|21[) is then 
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Fig. 9. A gluon ladder starting from a quark feo- The quasireal emissions have 
momenta qt, while the virtual links arc denoted k%. 



replaced by P : 



I(l + z 2 ) _ h i(l + z 2 ) , 3 



P = )) ' + -6(1 - z), (23) 



1-z 2 
where the definition of the "+ prescription" is given by 



f(z) =[ dz mzm (24) 



(i - z)+ w ; (1-3) 

for an arbitrary function f(z). We note in particular that this definition 
implies that 

r p(z)dz = 0, (25) 



which guarantees that the number of quarks is conserved. 

4-4- Exclusive final states 

We want to emphasise that eq. (|22|) only describes the probability for an 
interaction, that is the total photon-proton cross section expressed in terms 
of the structure function F2 ■ To find the properties of exclusive final states 
it is necessary to add final state radiation within angular ordered regions, 
as indicated in fig. |SJ3. Just as the final state emission in e + e~-annihilation, 
these emissions do not affect the total cross section, and should therefore 
also be associated with appropriate Sudakov form factors. 

4-5. Gluon rungs and small x 

Also gluons can emit gluon radiation. Study the emission of gluon qi in 
fig. With the notation in this figure the eikonal current corresponding 
to the first term in eq. (fT2*|) gives g{k\ + ^2)^2 ( e fci e fc 2 )- This term gives 
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a contribution similar to that for emission from a quark in eqs. (|171 118|) , 
Including now only the term singular at z = 1 we get 

3a s dzj dg 2 u 

-, ^%±,i - q±.i-i) ■ S. (26) 

7T 1 - Zi q± ti 

The spin factor efc 2 implies that the gluon spin is conserved along the 
vertical line in the ladder in fig. 03 We also note that, just as for quark 
emission, the pole at z = 1 implies that z-values close to 1 dominate. Thus 
gluon &2 can be regarded as identical to the recoiling parent k\\ it inherits 
both the spin and most of the energy, while the gluon q2 has to be regarded 
as a newborn soft gluon. 

The non-eikonal contribution to the current is proportional to the re- 
maining terms in eq. (JHJ), ~(fa + q 2 )e k2 {^k^q 2 ) ~ (<?2 - k 2 )e kl (efc 2 e q 2 )- Here 
the first term dominates (for a conventional gauge choice), and gives a con- 
tribution ~ dzi/zf. 

(27) 



7T 



Thus the "soft" gluon goes up the ladder, and it is the (quasireal) gluon 
q2, which takes over both the spin and most of the energy of the parent k\. 
The soft daughter does not replace the mother, and therefore there is no 
Sudakov form factor in leading order. 

We note that soft gluon links can also be emitted from quark legs, as also 
indicated in fig. |§J The probability for this process has a different colour 
factor, and is proportional to 

4a s dz dq\ 

l^T^f' (28) 

and also in this contribution there is no Sudakov form factor to leading 
order. 

For gluonic chains the sum of the expressions in eqs. (|26jl and (|27|) replace 
the corresponding factor in eq. (|2U|) . The non-eikonal terms in eq. (|27|) give 
the dominant contributions for small z-values, and therefore also for small 
x = (II Zi) xq. Thus gluon ladders are most important for the growth of the 
structure functions for small x, and to leading order in In 1/x and InQ 2 the 
gluon density can be written 

G(x, Q 2 ) - £ fl { / - q±,i-i)\ S(x-f[ z r x )9(Q*-qi 

n j iy Z i ?±,i J j 

(29) 
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The gluon ladder may start from an initial quark with the coupling in 
eq. (|28[). and in DIS it must also have a quark link at the end, as the photon 
only couples to the electrically charged quarks. It is, however, the properties 
of the dominating gluonic ladder, which determines the asymptotic growth 
rate for small x. 



4-6. Double Leading Log approximation 
Assume that both Q 2 and 1/x are very large. We introduce the notation 

f[ 2j (30) 



Pprot 



on 



a = (31) 

7T 

«i = HqljK 2 ) (32) 
k = ]n(l/xi) (33) 

For a fixed coupling a, we then find from eq. (|29j) 

g ~ e„. {n? i lnQ2 adKiOfa - • mi lal/x dw{u - u-i)} = 

_ V -n (InQ 2 )" (lnl/z)" _ 

= J (2^ InQ 2 In 1/x) ~ exp (2^/a lnQ 2 lnl/x) , (34) 



where Jo is a modified Bessel function. This result corresponds to the Double 
Leading Log (DLL) approximation. 

For a running a s we define the parameter «o by the relation 

a = H&m- (35) 

We then have instead of dq 2 _ i /qj_ i = dKi a factor dni/ni = dhiKj, which 
gives the result 



G~exp 2i/a • InhiQ 2 • lnl/x . (36) 



5. Small x, the BFKL region 

5.1. Non-ordered ladders 

Now assume that Q 2 is not large, while x is still kept small. In this case 
the gj_-ordered phase space is small. Therefore q_|_-non-ordered contributions 
are important, even if they are suppressed. 
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We mentioned after eq. (12011 that for a (/j_-ordered chain the transverse 
momentum of the (quasi)real emissions, qj_i, are approximately equal to 
the transverse momentum of the virtual links, denoted j in figs. |H] and 
Therefore the variables q±^ in eqs. (f2"U|) and (|2U)l could be replaced by 
the link variables k±^. This is no longer the case for non-ordered chains. 
Transverse momentum conservation implies that qj_ j = kj_ j_i — kj_ j, and 
therefore we have, provided kj^-i and kj_ j are not approximately equal, 

Q±,i ~ max(fei,i»*:i,i-i)- ( 37 ) 

Emissions for which kj_ j_i ~ kj_ j will be discussed in subsection 5.5. They 
give no contribution to the structure functions, and can be treated as final 
state radiation. For other emissions eq. (|^7jl implies that if e.g. is larger 
than the neighbouring links and fci^+i), then q±j « (Ll^+i ~ fcj^j. 

If on the other hand k± t i is smaller than its neighbours, then its value is not 
close to any of the transverse momenta q±. Thus the quasi-real momenta 
Q±,i are approximately determined by the virtual links k±i, but the reverse 
is not true (without the knowledge of the azimuthal angles of all the vectors 
q_L,i)- To specify the chain we therefore have to specify the link variables 
k± t i, and the distinction between q±^ and k±^ is essential. 

Classical bremsstrahlung due to a short impulse during a time At ~ 1 /Q 
contains the factor 



'i/Q 



dte 



iwt 



hju/Q - 1) 




for oj « 
%x, for oj » 



(38) 



Thus for oj » Q there is a relative suppression by a factor Q 2 /oj 2 . In a 
relativistic calculation we should make the replacement oj — > q±. We then 
find that for q 2 ^_ > Q 2 the emission is not totally excluded; it is only reduced 
by a suppression factor Q 2 jq\- 

Assume that we have an ordered chain up to the last link with k\ ~ 
q\ n , which can be smaller or larger than Q 2 . We then find for these two 
cases: 



k ±,n > k ±.n-V and n ° 



kl, n < Q 2 - 

In this ordered case we have Q 2 > q\ n 
extra suppression. According to eqs. 1)20(1 and (|2"§)l the contribution 
to the structure function from these chains contain a factor l/qj_ n - 
Therefore we have 

a ~7v- F ~7W^- ( 39 ) 

Q Q T±,n 
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(a) (b) (c) 

Fig. 10. (a) A situation where k\ n is larger than the virtuality Q 2 of the probe 
can be interpreted as a boson-gluon-fusion event, i.e. a hard subcollision between 
the probe and gluon fc n -i- (b) If k\ n _ 1 is the largest virtuality, we have a hard 
subcollision between gluons k n _ 2 and k n . (c) When the gluon with the largest 
transverse momentum is located in the middle of the chain, it corresponds to a 
hard subcollision between a proton and a resolved photon. 



• > Q 2 - 

This implies Q 2 < q\ n « k\ n > k\ n _ v and in this situation there is 
a suppression factor Q 2 /q\ n - Thus we have 

i i Q 2 i 

9±,n «_L,n K ±,n 

In the latter case we see that the process can be interpreted as a hard 
subcollision between the probe and the gluon k n _i, cf. fig. ITUb. (Instead of 
a photon we imagine here a hypothetical colour neutral probe, which can 
interact directly with gluons.) We recognize the expected result from a hard 
scattering with momentum exchange t = —k\ n , which is then the largest 
virtuality in the process. 

Assume now that we have a situation where the last but one gluon link, 
k n -\, has the largest transverse momentum, which implies: 

• Q 2 < n < q\^i ~ fei >n _l ~ q±, n -2 > ^±,n-2- 

We expect here a suppression factor Q 2 /q± n _i ~ ''/k 2 _ n _ 1 . This fac- 
tor can be written as a product of two separate factors, one for each step 
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downwards in k±: 



« 2 - *~ * (41) 



h 2 h 2 h 2 



As illustrated in fig. HUb . this process can be interpreted as a hard sub- 
collision between the links k n (now coming from above) and k n -2, with 
momentum transfer k\ 1 . The outgoing gluons have approximately equal 
but opposite transverse momenta, g_i_ in _i and q± jH , and the cross section 
satisfies the expected relation a ~ 1/^1 n-l- 

These results can be generalized to the situation in fig. 110b . where k± 
increases continously from the proton, up to a maximum value fc_L )max > an d 
then decreases in k± down to the virtuality, Q 2 , of the photon. This pro- 
cess has a weight with a factor for every links except /c max , which 
instead gives a factor l/k± max - We note here that, although upwards and 
downwards steps are treated differently, the final result is symmetric in the 
sense that we could equally well have started the chain from the photon 
end, and proceeded towards the proton. Therefore it is identical to the re- 
sult obtained from a DGLAP evolution for both the proton in one end and 
a resolved photon in the other end, up to a central hard scattering, with 
momentum transfer /cj_ jmax . 

With enough energy it is apparently also possible to have ladders in 
which the transverse momentum goes up and down with two or more local 
maxima. When expressed in the virtual links k±, each step downwards 
corresponds to a suppression factor k'j_ i /k 2 _ i _ 1 . The net result is a factor 
l//cj_ for every local maximum k±, but no fcj_-power for links corresponding 
to a local minimum. 

We formulated the result in terms of the link variables k\ i , because 
these could be interpreted as independent variables. Transverse momentum 
conservation implies that the quantities q± i are not independent, but con- 
strained by the relation q\ max « q\ max+1 , and they do not fix the value of 
k\ for a link which represents a local minimum. We note, however, that the 
weight for the chain corresponds exactly to the product Ylf q±\- The factor 

^1 max ' 1± max+i & ves the factor l/kj_ max , and the k± for a local minimum, 
which does not equal any q±, does not appear in the weight. 

We further note that, if the azimuthal angles are not averaged over, but 
properly accounted for, we have d 2 k± : i = d 2 q± : i. As a consequence we see 
that our result exactly corresponds to the ordered result in eq. (|29|). if we 
make the replacement 

d Q±,i d 2 q ±yi 



omit the ordering ^-functions and instead include a factor 



(42) 
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9{q±,i — m i n (^_i_,i; &±,i-i))- This factor excludes soft emissions for which 
k_i_,i ~ k^^-i, as mentioned above and is further discussed in section 5.5. 
Thus the expression 

G ~ i / ^- — d q ^ l 9(q± : i - mm(k ±ii , fcjL,i-i))> S(x - \\zj x ) 

n i V d7F Z * ^±,1 J j 

(43) 

gives a proper description, both for large Q 2 and ordered chains, and for lim- 
ited Q 2 when non-ordered chains are important. To get the non-integrated 
distribution function Q, we just have to add a factor 5(k± tTl + X)? Qi-j)> 
which follows from conservation of transverse momentum. 

This result also demonstrates the symmetry discussed above. An im- 
portant consequence of this symmetry is that the formalism also is suitable 
for describing hard subcollisions in hadronic collisions. This will be further 
discussed in section 6. 



5.2. Effective phase space 

With more than one local maximum in the k± chain, we find for ev- 
ery step down from k'j_ i _ 1 to kj_ i a suppression factor k\ i /ki Li _ 1 . Ex- 
pressed in the logarithmic variable k = lnkj_, this corresponds to a factor 
exp[— (Kj_i — Ki)}. This implies that the effective range allowed for down- 
ward steps corresponds to approximately one unit in k. Consequently we 
find instead of the phase space limits in eq. (JT5|) , the following boundaries 
(replacing q'j_ i by k\ i} which are the relevant variables for non-ordered 
chains) 

lnk 2 ±ji > lnfci^-l. (44) 

This modification changes the DLL result in eqs. (|34j) and 1)36(1 in a quali- 
tative way, as described in the following subsections. 



5.3. Fixed a s 

For a fixed a s , and writing k for hiQ 2 , we then find for the transverse 
momentum integrals in eq. (|34[) instead of 

/ ldK i 0(K i -K i - 1 ) = ?-, (k = InQ 2 ) (45) 
Jo n! 

the following result 

/ n^c^-^-i- 1 ) ^ , ( 46 ) 
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When k is very large we recover the DLL result in eq. (|H4J) . but for 
smaller values of k we find instead using Sterling's formula 



which implies 



k small ~ ^ ~ e n (47) 

n! n! 



G ^ [aeln(l/x)] w = ^ a ln(1/x) = J_ 
^ n! x A 



with 

A = ea«2.72a. (49) 

To obtain the non-integrated distribution Q(x,k'j_), we add the factor 
<5(kj_ — kj_ )n ), but this does not modify the qualitative result, if k now 
denotes lnAr[. 

The result in eq. (|48j) is relevant for ln/c^ < ealn(l/x). In this range 
the chain corresponds to a random walk in ln/c^. Our result should be 
compared with the result from the leading order BFKL equation, which 
gives 

A = 41n2a « 2.77a. (50) 

We see that this simple semiclassical picture describes the essential features 
of BFKL evolution. 



5.4- Running a s 

For a running a s the steps in transverse momentum are determined by 
the factors 

a o-j — To — ~ = cxodCln In k\ A (51) 
lnk ±,i 

When k± is large, one extra unit in In k\ corresponds to a very small extra 
space in lnlnfc^. This implies that, once the chain has reached large k±- 
values, it is very difficult to come down again. Therefore dominant chains 
will contain an initial part with low kj_ and steps up and down, and a second 
(DGLAP-like) part with increasing k± up to the final k±- and x-values. As 
a consequence the structure functions can be well described by DGLAP 
evolution from adjusted input distributions fo(x, &j_o)' which grow for small 
values of x. 

For further discussions of the results in this and the previous sections, 
see ref. [XT] , 
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5. 5. Exclusive final states 



As discussed above in connection with large-Q 2 phenomena, final state 
radiation must be added, in order to describe the properties of exclusive final 
states. The kinematic regions in which such emission should be allowed is in 
general dependent on the formalism used. The quasireal gluons, denoted qi 
in figs. |H1 and EJ are allowed to emit final state radiation in angular ordered 
regions, but besides within these regions further soft emissions may also be 
possible. 

In all discussions of the results in the previous sections we studied chains 
in which at every vertex the emitted q± approximately equals the larger of 
the adjoining k±. We could also imagine the emission of a softer gluon, 
with qx,i « k± i—i, which then implies kj_ : j ~ kj^-i. From the discus- 
sion about formation time in section 2.2, such emissions should take place 
at a longer time scale, after formation of the interaction chain in figure |SK 
or El They should therefore not influence the reaction probability, i.e. the 
structure function, but only affect the properties of exclusive final states. 
Consequently they may be treated as final state radiation, and must be as- 
sociated with appropriate Sudakov form factors. This is indeed the result 
obtained in a perturbative calculation, in which such emissions are com- 
pensated by virtual corrections, and in the BFKL formalism this effect is 
described by treating the link gluons as Reggeized gluons. Such final state 
radiation may be emitted both from gluon and quark links, and one example 
was indicated in fig. |Hj3. 

These emissions are also treated as final state radiation in the Linked 
Dipole Chain (LDC) model for DIS [I], which will be further discussed in 
the next section. We note, however, that the separation between initial and 
final state radiation is not defined by nature, but depends on the formalism 
used. This separation is defined in a different way in the CCFM formalism 
P], in which some emissions, for which kj_ j ~ kj_ are treated as ini- 
tial state radiation. This implies that in this formalism the initial chains 
generally contain a larger number of links, and that final state radiation 
is correspondingly allowed in a more restricted kinematic region. A conse- 
quence is that a larger set of chains contribute in the calculation of Q, with 
each chain given a smaller weight. This reduction is represented by a "non- 
eikonal form factor" , in such a way that the BFKL result is reproduced to 
leading log accuracy. For further discussions of this difference between the 
formalisms see refs. [I] and |12j . 
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6. Linked Dipole Chain model 

6.1. LDCfor DIS 

The Linked Dipole Chain (LDC) model is a reformulation and gen- 
eralization of the CCFM model. Thus it is based on perturbative QCD 
calculations, and can be regarded as a formal quantum mechanical deriva- 
tion of the semiclassical results presented above. 

The LDC model is based on the observation that the dominant features 
of the parton evolution is determined by a subset of emitted gluons, which 
are ordered in both positive and negative light-cone components, and also 
satisfy the relation 

q±j > min(/c_ Lii , fc_L,j_i). (52) 

In LDC this subset of "primary" gluons forms a chain of initial state radi- 
ation, and all other emissions are treated as final state radiation (including 
those emissions, which would be emitted from a virtual link ki, rather than 
from a quasireal emission % as discussed in section 5.5). In ref. [I] it is 
shown that adding the contributions from all different CCFM-chains with 
the same primary gluons, with their non-eikonal form factors, gives the 
following simple result: 

where the link momenta kj^j = *LL,j are determined by transverse mo- 
mentum conservation. This implies that 

i±,i ~ max (^i,i) ( 54 ) 

From this relation we find for a step up or down in k±: 

• Step up in k±: k±^ > k± t i-i =4> q± t i ~ k±^, which implies 

d 2 q±,i d 2 k\i 



• Step down in k±: k± t i < fei^-i => q±A ~ which implies 

d 2 q±,j _ d 2 k ±ti _ d 2 k ±ti k\- 

„2 ~ L.2 L.2 h.2 

1±,i K XA-1 K ±A ^±,1-1 



(55) 



(56) 



We here recognize the suppression factor k\ i jk\ i _ x in eq. (|1T|) . associ- 
ated with a step down in k±. 
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Fig. 11. The description of F 2 data as a function of x for different Q 2 . To separate 
the results, what is shown is F 2 + i for Qf . Both small and large x data from 
HI, ZEUS, NMC and E665 are included QSj. The favoured fit is the one called 
standard and denoted by dotted lines. In this fit both gluon and quark links are 
included. For a detailed description of the fits, see ref. |15| . 

Non- leading contributions from 1/(1 — z)-poles and quark links are added 
with Sudakov form factors. In this formalism it is natural to include a 
running coupling 

= a s (max(fci i , fci^x)) (57) 

Note, however, that this increase of a s for soft emissions necessarily 
implies that a cutoff is needed for soft k± [El- This dividing line between 
the perturbative and the nonperturbative regimes has to be adjusted by fits 
to experimental data. 

A Monte Carlo implementation, called LDCMC, is developed by H. Khar- 
raziha and L. Lonnblad |T3|. This program does reproduce a large set of 
experimental data. Fig. ^2 shows a fit to F2 ^2] , compared with data from 
HERA and fixed target experiments. The corresponding gluon distribution 
is shown in fig. 1121 and we note that this result shows good agreement with 
the results from the CTEQ [T7] and MRST Q2] fits. (Note that the LDC re- 
sult is fitted to DIS data only, while the CTEQ and MRST fits include data 
from hadronic collisions.) As mentioned above, the LDC model is a gener- 
alization of the CCFM model, which is implemented in the event generator 
Cascade [H]. Both ldcmc and Cascade reproduce the main features of 
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Fig. 12. The LDC integrated gluon distribution function for two different fits (called 
"gluonic" and "gluonic-2" ) , compared to the corresponding results from Cascade 
denoted JS (dash-dotted curve), CTEQ (short-dashed curve) and MRST (long- 
dashed curve), for (a) Q 2 = 16 GcV 2 and (b) Q 2 = 100 GcV 2 . 



the final state properties [20]. A problem is, however, that both models are 
able to reproduce the forward jet cross section only if the non-singular terms 
in the splitting functions are omitted [20]. We do not know whether these 
terms are compensated by some dynamical mechanism, or if the modelling 
of the proton fragmentation has to be improved. 

6.2. Hadronic collisions 

The LDC formalism can be applied also to hadron-hadron collisions 
\'21\ 122 j . In hadronic collisions multiple interactions are needed to describe 
associated multiplicity and transverse energy in events with high p± jets. 
Such multiple interactions can originate either from two hard subcollisions 
in a single parton chain, or from more than one chain in a single event. 

The LDC formalism is very suitable for describing such events. A po- 
tential problem arises because good fits to DIS data can be obtained for 
different values for the low k± cutoff, provided that the soft input structure 
function is adjusted accordingly. This is important because the number of 
hard chains depends on this soft cutoff. If the cutoff is increased a more 
singular input gluon distribution is needed in the fit to 7*2- These soft input 
gluons, present at scale Qq, may also form "soft chains" between the two 
colliding protons. These soft chains then have no emitted gluons with q± 
above Qq, but they do produce hadrons with low p±. From the fits to DIS 
data we then find that the reduction in the number of hard chains is just 
compensated by an increase in the number of soft chains, as is illustrated 
in fig 1131 As a result we see that, in this formalism, the (average) total 
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Fig. 13. The cross section per chain in the LDC model as a function of the cutoff, 
k\ n . The dashed line is the cross section for chains with at least on emission above 
the cutoff, the dotted line is for soft chains without emissions above the cutoff, and 
the full line is the sum of the two. Note that the input parton densities have been 
re- fitted for each value of k±o- 

number of chains in a pp collision can be fixed by fits to DIS data. This re- 
sult apparently implies a very strong connection between DIS and hadronic 
collisions |22j . 

This formalism should also be applicable in studies of effects of unitari- 
sation and saturation. Separate branches of a chain can interact with a 
colliding hadron, or different chains can join each other, and these possibili- 
ties are presently under study. Also applications to diffractive collisions are 
of interest. 

7. Conclusions 

We have seen that the behaviour of DIS in the DGLAP and BFKL 
domains can be understood in a semiclassical intuitive picture. A more 
quantitative description is offered by the LDC formalism, which smoothly 
interpolates between these two kinematical regions. 

The formalism also offers a link between DIS and hadronic collisions. 
Within the formalism fits to data on F2 in DIS can be used to predict the 
number of chains and hard subcollisions in a high energy proton-proton 
collision. 
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Further work on multiple interaction and saturation is in progress, as 
well as studies of diffractive scattering. 
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